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Abstract 

The control of the spread of dengue fever by introduction of the intracellular parasitic bacterium 
Wolbachia in populations of the vector Aedes aegypti , is presently one of the most promising tools 
for eliminating dengue, in the absence of an efficient vaccine. The success of this operation requires 
locally careful planning to determine the adequate number of individuals carrying the Wolbachia 
parasite that need to be introduced into the natural population. The introduced mosquitoes are 
expected to eventually replace the Wolbachia-hee population and guarantee permanent protection 
against the transmission of dengue to human. 

In this study, we propose and analyze a model describing the fundamental aspects of the com¬ 
petition between mosquitoes carrying Wolbachia and mosquitoes free of the parasite. We then use 
feedback control techniques to devise an introduction protocol which is proved to guarantee that the 
population converges to a stable equilibrium where the totality of mosquitoes carry Wolbachia. 
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1 Introduction 

1.1 Arboviroses and vector control 


Arboviruses (arthropod borne viruses) are viruses transmitted to humans by arthropod, such as the 
mosquito. They are pathogens of many and important diseases, putting at risk considerable portions of 
the human population, and infecting millions of people every year. Mosquitoes ( Culicidae family of the 
insects) are a huge public health concern as they are vectors of many arboviroses such as yellow fever, 
dengue and chikungunya. 

The control of these diseases can be achieved by acting on the population of mosquitoes, and in absence 
of vaccine or curative treatment, it is essentially the only feasible way. Application of insecticides for both 
adults and larvae and mechanical removal of breeding sites are the most popular strategies to control the 
population of mosquitoes. The intensive use of insecticides, however, has negative impacts for humans, 
animals and the environment. Besides, the gradual increase of mosquito resistance to insecticides usually 


leads to partial or complete decrease of the efficiency of this strategy in the long run (de Freitas and 


Valle 2014 Montella et al 2007). In addition to chemical control and mechanical removal of the breeding 


sites, alternative or supplementary vector control strategies have been proposed and implemented, such 
as the release of transgenic or sterile mosquitoes (Alphey et al 20l0 Alphey 2014). Notice that an 


intrinsic weakness of the techniques listed above lies paradoxically in the fact that they aim at the local 
eradication of the vector, whose disappearance offers no protection against subsequent reinvasions. 

Recently the release of Aedes aegypti mosquitoes infected by the bacterium Wolbachia has been 
proposed as a promising strategy to control dengue and chikungunya, due to the fact that this bacterium 
limits severely the vectorial competence of Aedes aegypti. The international program Eliminate Dengue 
(Hoffmann et al) is currently testing in the field this strategy, in several locations around the world: 
Australia, Indonesia, Vietnam, Colombia and Brazil. The release of infected mosquitoes with Wolbachia 
has the advantages of being safe for humans and the environment and inexpensive when compared to 

M0). 


other control strategies (Popovici et al 


1.2 Wolbachia sp. as a biological control tool 


Wolbachia sp. is a genus of bacteria that is a common intracellular parasite of many species of arthropods. 
It is often found in anthropophilic mosquitoes such as Aedes albopictus or Culex quinquefasciatus but 


there is no report of Aedes aegypti naturally infected by this bacterium (Rasgon and Scott 2004). 


There is evidence that the spread of certain strains of Wolbachia in populations of Aedes aegypti 


drastically reduces the vector competence of the mosquito for dengue and other diseases (Moreira et al 


2009 Blagrove et al 2012 2013). Some strains of Wolbachia reduce the lifespan of the mosquito, 


consequently limiting the proportion of surviving mosquitoes at the completion of the incubation period. 
More importantly, Wolbachia appears to decrease the virulence of the dengue infection in the mosquitoes, 
increasing the incubation period or blocking the virus, which also reduces the overall vector competence. 

The infestation of natural Aedes aegypti populations by Wolbachia-cont&mm&ted strains can be 
achieved by releasing in the field a large number of Wolbachia-iniected mosquitoes bred in laboratory. 


Experiments have been conducted successfully in Australia (Hoffmann et al 2011), Vietnam, Indonesia 


and are currently being applied in Brazil, within the Eliminate Dengue Program. In these experiments, 
the introduction of a number of Wolbachia- infected mosquitoes in the population triggered a Wolbachia 
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Uninfected cf 

Infected cf 

Uninfected 9 

Uninfected 

Sterile eggs 

Infected 9 

Infected 

Infected 


Table 1: Cytoplasmic incompatibility and vertical transmission of Wolbachia bacteria. The state of the 
offspring is indicated, depending on the parents status 


outbreak whose outcome was the fixation of the bacteria in the population, with more than 90% of preva¬ 
lence. The effects of this fixation of Wolbachia on the dynamics of dengue in the field is currently under 


investigation, but preliminary results are encouraging (Frentiu et al 2014). If tractable in practice, this 


method has certainly the advantage of offering certain resilience to subsequent invasion of Wolbachia-iree 
mosquitoes. 

Several mathematical models of the dynamics of invasion of Wolbachia in a population of mosquitoes 


have been proposed, with distinct objectives. For example, Turelli (2010) describes a simple model with 


a single differential equation, sufficient to reveal the bistable nature of the Wolbachia dynamics. Models 


for spatial dispersion are analyzed 

in (Barton and Turelli 

2011 

) and (Hancock and Godfray 

2012 

). In 

(Hughes and Britton 

2013 

Ndiietal 2015 

), models are presented that assess the effect of the Wolbachia 

in dengue dynamics. 

Koiller et al 

(2014 

) describes a data driven model suitable to estimate accurately 


some biological parameters by fitting the model with held and lab data. The model to be presented here 
is a simplified version of the latter. 


1.3 Description of the problem 

A key question about the introduction of Wolbachia in wild mosquitoes concerns the effective strategies 
of release of infected mosquitoes in the held that can be applied with limited cost to reach the desired 
state of 100% of Wolbachiarcarrying mosquitoes. In this paper we propose and analyze a simple model 
of the dynamics of Wolbachia , that allows to investigate these strategies. 

The main features of the natural dynamics of Wolbachia that have to be present in the model are the 
vertical transmission and the peculiar interference on the reproductive outcomes induced by cytoplasmic 
incompatibility (O’Neill et al 1998). The transmission of Wolbachia occurs only vertically (i.e. from 


mother to the offspring), there is no transmission by contact. Cytoplasmic incompatibility (Cl) occurs 
when a female uninfected by Wolbachia is inseminated by an infected male, a mating that leads to 
sterile eggs. It provides a beneht to infected females against uninfected ones and therefore facilitates 
the Wolbachia spread. The strains wMel and wMelPop of Wolbachia that are being used in the field 
experiments with Aedes aegypti induce almost total Cl (Walker et al 2011). Table [l] schematizes the 
results of the mating of infected and uninfected mosquitoes when the Cl is 100%. The model that we 
propose below captures all of these features and is simple enough to allow for a fairly complete analysis. 

The corresponding system is shown to possess two unstable equilibria, which correspond to extinction 
of the two populations and a coexistence equilibrium; and two locally asymptotically stable equilibria, 
which correspond to Wolbachia-free and complete infestation equilibria. The release strategy we propose 
here is based on techniques from Control theory — the released quantity of infected mosquitoes being 
the control input. Using the fact that measurements are achieved and available during the whole release 
process, we propose a simple feedback control law that uses this knowledge to compute the input value. 
The key result of the paper (Theorem [l2|) proves that this control law has the capacity to asymptotically 
settle the bacterium from whatever initial conditions , and in particular from the completely Wolbachia- 
free equilibrium. A major advantage of feedback compared to open-loop approaches (when the release 
schedule is computed once for all before the beginning of the experiment), is its ability to cope with 
the uncertainties on the model dynamics (e.g. in the modeling of the life stages and the population 
structure), on the parameters (population size, mortality, reproductive rates, etc.), and on the conditions 
of the realization (in particular on the size of the population to be treated). 

Up to our knowledge, the present paper constitutes the first attempt to use feedback approach for 
introduction of Wolbachia within a population of arthropods. Notice that we treat here only the case 
of the release of Wolbachia-positive larvae and full information on the quantity of Wolbachia-negative 
larvae. Yet, the same dynamical model offers the ability to study other configurations, both for the 
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Notation 

Meaning 

ay, aw 

V 

H 

/ik 

Hu, Hw 

Fecundity rates of uninfected and infected insects 

Rate of transfer from the preliminary to the adult stage 

Mortality rate of uninfected and infected insects in preliminary stage 
Characteristic of the additional mortality rate in preliminary stage 
Mortality rates of uninfected and infected insects at adult stage 


Table 2: List of parameters of model 0 


control and the observation, and the corresponding issues will be examined in future work. 

The paper is organized as follows. The simple model used in the sequel is introduced in Section [2j 
and normalized. The analysis of the uncontrolled model is made in Section [3j showing the bistability 
announced above between Wolbachia-iree equilibrium and full infestation. A proportional control law 
is then proposed in Section [4j and proved to lead to global stability of the full infestation equilibrium. 
Simulations are provided in Section [5] Last, concluding remarks complete the text in Section [6] 

Notation For n € N, we let R 71 denote the n-dimensional Euclidean real space, and with R" we refer to 
the cone consisting of vectors in R™ with nonnegative components. We write max{a; b} (resp. min{a; b}) 
for the maximum (resp. minimum) of two real numbers a,b. 


2 A simple model of infestation by Wolbachia 


The simplified compartment model we introduce includes two life stages: a preliminary one, gathering 
the aquatic phases (egg and larva) where the mosquitoes are subject to space and food competition; 
and an adult one, representing all the posterior aerial phases (pupae, immature and mature adult). 
Accordingly, we will denote L and A the corresponding state variables. The uninfected and infected 
(by Wolbachia ) populations will be distinguished by indexes U and W respectively, so we end up with 
a four state variables model, namely Ly,Lw and Ay, Aw, that represent the numbers of uninfected, 
resp. infected, vectors in preliminary and adult phases. 

We propose the following evolution model. 

Ay 

Ly = Qy--—-Ay — foLy — fl( 1 + fc(Lyv + Ly))Ly (la) 

Ay + A w 

Ay = vLjj — /iyAy (lb) 

Lw = CX-W^-W ~ vLw — m(1 + k(Li\y + Ly))Lw + u (lc) 

Aw = i'Lw — Hw-^-w (Id) 


All the parameters are positive, their meaning is summarized in Table [2j 

Most aspects of this compartment model are rather classical, we now comment on the most original 
modeling choices. The differences between the behaviors of the two populations lie in the different 
fecundity and mortality rates. The mortality during the larva stage and the duration of the latter are 
considered unmodified by the disease. Also, the (quadratic) competition term is assumed to act equally 
on both populations, with an effect proportional to the global number of larvae. 

The first effect of Wolbachia is to diminish fertility and life duration, leading to reduced fitness for 
the infected mosquitoes. This assumption will correspond to the choice of parameters made in ([5| below. 
The second effect of Wolbachia , namely the cytoplasmic incompatibility, modeled here as complete, 
forbids fecund mating between infected males and uninfected females. This is rendered in ( |la| ) by a 
recruitment term proportional altogether to the number and to the ratio of uninfected adults. Notice in 
this respect that the model apparently does not make distinction between males and females. In fact, 
one could introduce in place of the two variables Ay, Aw, four variables Fy, Fw, My, Mw representing 
the quantities of female and male adults. The recruitment terms in resp. fla| ), would then naturally 
be replaced by expressions proportional to Fw, resp. MAyfenyFy. However it is easy to see that the 
proportion between males and females remains constant in this more complex system, as long as the 
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mortality rates for healthy and infected insects are equal for the males and the females, and the passage 
from larvae to adults occurs with a constant sex ratio. One can therefore use a single variable to take 
account of the uninfected adults, and a single one to take account of the infected adults, just as done 
in 0- In other terms, provided the sex ratio is constant and the mortality is sex-independent, a sexual 
model yields no more information than 0 „ 

As a last comment, notice the term u in equation ( |Tc| ) . The latter is an input variable , modeling the 
on-purpose introduction of infected larvae in the system in order to settle Wolbachia. 


2.1 Normalization and general assumption 

In order to reduce the number of parameters and to exhibit meaningful quantities, we now normalize 
model Q. Defining 




_kp_ ] 

v + n 


t 


v + 
u{t) := 


A v (t ) := ~A V 


v + fj, 


r) = U,W 


kp 


(v + nY 

with the following choice of dimensionless parameters 

Mt) nyV _ 


t 


v + n 


'■= 


°' (i' + MK 


T) = U,W 


IS ~\~ f-l \is ~r 

the following normalized model is deduced, that will be used in the remainder of the paper. 

Au 


Lu = 


0 Atj + A\ 


-Au — (1 + Lw + Ljj)Lu 
i £7 - 1 - ri w 

Au = Lu - 'yuAu 

L\V = Aw — (1 + Lw + Lu)L\y + u 

Aw = Lw — jwAw 

The state variable for system @ will be denoted 

x := (Lu, Au, L w , A w ) , 
and for sake of simplicity, we write (J3| as 

x = f(x) + Bu , 

where / and B are defined as 


/ 0 ) : = 


/ luTtf Au + Aw A v (1 + L W + Lu)L v ^ 


(°\ 

Lu - luAu 

B — 

0 

IwLYq A w - (1 + L w + L L j)L w 

, D . 

1 

\ L w - 7 wAw / 


w 


We assume in all the sequel 


> n™ > 1 . 


( 2 ) 


(3a) 

(3b) 

(3c) 

(3d) 


(4) 


(5) 

Assumption ([5]), which will be valid in the remainder of the paper, ensures the sustainability of each of 
the two isolated populations, with an even better sustainability for the non-infected one. See Theorem 
[7] below for more details. 


3 Analysis of the uncontrolled system 

The uncontrolled system is obtained by taking zero input u , that is: 


x = /( x) 


( 6 ) 
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3.1 Well-posedness, positivity and boundedness 

One first shows the well-posedness of the Cauchy problem related to equation (Jh]) for nonnegative initial 
conditions. 

Theorem 1. For any initial value in there exists a unique solution to the initial value problem 
associated to system (§. The latter is defined on [0,+ 00 ), depends continuously on the initial conditions 
and takes on values in • Moreover, it is uniformly ultimately bounded. 

Remark 2. The previous result shows that system ([6| is positive. Therefore, when talking about 
“trajectories”, we will always mean trajectories with initial values in The same shortcut will be 
used for all positive systems considered later. 

Let us introduce the following definition of an order induced by a cone, that will be instrumental in 
proving Theorem [I] 

Definition 3. Let n £ N and K, C R n be a closed convex cone with nonempty interior. We use >ic to 
denote the order induced by K, that is: \/x,x' £ 

x >k. x l <t=> x — x' £ K. . 


Similarly, one defines 


x >/c x 1 <t=> x — x l £ K. and x A x l 
x x’ <t=> x — x' £ int K. 

As usual, one will write x <n x' to mean x' >ic x. 

Proof of Theorem^ 7} Function / in © is clearly well-defined and continuous in R+, except in points 
where Au = Aw = 0. Due to the fact that 0 < Al ^ff Aw < 1, the quantity Al A + Aw A\j tends towards 
zero when one approaches such points, and /( x) can thus be defined by continuity when Au = Aw = 0. 
In addition, the right-hand side is clearly locally Lipschitz in K^_, and classical result ensures the local 
well-posedness of the initial value problem, as long as the trajectory does not leave this set. 

The invariance property of the set is verified due to the fact that 

Va; £ R^_, V* £ {1,2,3,4} : aq = 0 => fi{x) > 0 . 

Let us now show that, for any initial condition in K^_, the associated trajectory remains bounded for 
all t > 0. With this aim, let us define 


L := Lu + Lw, 

7 := minlyc/jy^} > 0, 

Notice that, in view of hypothesis ([5]), 


A := Au + Aw, 

--L I U ’ ^0 

/•Co •— - 

7 


max{7,7^; 7 w^q } 


Tin > 1 . 


It turns out that 

L = Julio-. — 777 — A L t + j w li^ A w - (1 + L w + Lu)(Lu + L w ) 

Au + Aw 

< (luTio Au + jwTio A w ) - (1 + L)L 

< jlioA — (1 + L)L 


and 

Now, the auxiliary system 


A < L — 7 A 


L' = jli 0 A' - (1 + L')L', A' = L' - 7 A' 


(8a) 

(8b) 

(9) 


( 11 ) 
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(1988)) for the canonic order induced by the cone R/. One may 


is evidently cooperative (see Hirsch 

thus use Kamke’s theorem, see e.g. (CoppeI|, 1965] Theorem 10, p. 29) or Smith (1995), and compare the 
solutions of ([6]) (with L and A defined by |8j)) and ©• One deduces 

L(t) < A(t) < A'(t), for all t> 0, 

whenever the solutions are considered with the same initial conditions. 


It may be shown without difficulty that system (11) possesses exactly two equilibria, namely 


1 


Z* :=(!,*, A*) := (0 ,0) and x* := (L*, A*) := ^H 0 - 1 ,-ft 

Due to linearization around each point shows that ir* is locally unstable, while x* is locally asymp¬ 
totically stable (LAS). On the other hand, notice that x * < R ^ x*. Using the local stability information, 
application of (Hirsch, 1988', Theorem 10.3) then shows that the stability of x* is global in the topological 


interior of R]_, and that this point is in fact attractive for any initial point distinct from x * = (0,0). In 
any case, all solutions of © converge to the order interval 

[**; z*] R 2 + := {x' 6R 2 + : L* < L' < L*, A* < A' < A*} . 

Coming back to the solutions of Q, the comparison method mentioned above now allows to deduce 
the same property for L , A defined in Q. Using finally the fact, proved in Theorem[lJ that the trajectories 
remain in R]_, the same bounds apply componentwise to Ljj, Lw and Ay, Aw respectively. In particular, 
all solutions of (J6| are uniformly ultimately bounded, and this yields global existence of solutions, and 
hence the proof of Theorem [l] □ 


3.2 Monotonicity 

One shows here that the uncontrolled system ([6]) is monotone. For sake of completeness, we recall here 
the definition of monotone and strongly order-preserving semiflows defined on a topological space X 
partially ordered by an order relation <£ generated by a cone K, (see jSmith (1995)). 

Definition 4. The semiflow $ is called monotone if 


$t{x) <K (x 1 ) whenever x </c x' and t > 0 . 

The semiflow 4 ’ is called strongly order-preserving i/4> is monotone and, whenever x </c x', there exist 
open subsets D, O' of X with x £ Q, x' £ Cl', and t > 0 such that 


< K $*( 0 ') , 


this meaning z </c z! , for all z G 4‘t(D), z' £ $t(t2 / ). The semiflow 4> is called strongly monotone if 4> 
is monotone and 

4>t(a;) <Ck ^(a/) whenever x </c x' and t > 0 . 

A dynamical system is said to have one of the properties above if its associated semiflow does. 

We now examine system § at the light of these properties. 

Theorem 5 . System is strongly order-preserving in R]_ for the order induced by the cone 

t:=Lxl_x R + x M + , (12) 


(that is such that: x x' <t=> Xi < x[, i = 1,2 and Xi > x[, i = 3 , 4 / 

System © is therefore monotone in R]_, but not strongly monotone, due to the fact that the trajec¬ 
tories departing inside the sets {x £ : Ly = 0, A v = 0} and {x £ R/ : L w = 0, A w = 0} remain 

in these sets and, consequently, do not verify strict ordering property for the two null components. 

Before proving Theorem [5j we summarize in the following result the behavior of the trajectories in 
relation with some parts of the boundaries. 
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Lemma 6. Let xq £ R{_. Then exactly one of the four following properties is verified by the trajectories 
departing from xq at t = 0. 


* x = £ 0,0 (that is, x(t) = £ 0 , 0 ; Vi > 0). 

* Aw = 0 and Ajj(t) > 0, Vi > 0. 

* Au = 0 and Aw{t) > 0, Vi > 0. 

* Aw(t) > 0 and Ajj(t) > 0, Vi > 0. 


Proof of Lemma^Q Clearly, one sees from (3b) (resp. ( |3d| ) that Au = 0 (resp. Aw = 0) if and 
-A[/(0) = 0 and Lu = 0 (resp. Aw( 0) = 0 and Lw = 0). Therefore, if (Au(0),Lu(0)) ^ (0,0) 
(Aw(0), Lw{0)) ^ (0,0)), then Au(t) > 0 (resp. Aw it) > 0) for all t > 0. This proves Lemma [6] 


only if 
(resp. 
□ 


Proof of Theorem [5| We now introduce the gradient of /. At each point x = (Lu , Au, Lw, Aw) £ R+ 
such that Au + A w > 0, V/(x) is equal to 


<-\-2Lu-L w 7^(l ~ {Au A h wY ) 

1 -7 u 

—Lw 0 

\ 0 0 


-Lu 

0 

— 1 — Lu — 2 Lw 
1 


\ Au t U A w )^ 

0 

7 wK 
~lw ) 


(13) 


Notice that, as a corollary of Lemma |6j either x = £ 0 ,o> or Au(t) + A w (t) > 0 for all t > 0. Therefore, 
the gradient can be computed at any point of a trajectory, except if the latter is reduced to £ 0 , 0 ■ 

For any x £ R^_, one verifies easily that 


V(i,j) £ {1,2} x {3,4}, ||(£) < 0, ^L(x) < 0, 

V(i,i) £ {1, 2} 2 U {3,4} 2 , i^j, Jg(£)>0. 


Hence the system is monotone. 

Moreover, except when Au = 0 or Aw = 0, the Jacobian matrix in (13) is irreducible, and the 
semiflow related to system § is therefore strongly monotone therein. On the other hand, trajectories 
confined to one of the sets {x £ : Lu = 0, Au = 0} and {£ £ : Lw = 0 ,Aw = 0}, also 

verify strong monotonicity, for the order relation restricted to the two non-identically zero components. 
These two remarks, together with Lemma [6j show that overall the strongly order-preserving property is 
verified. This completes the demonstration of Theorem [5] □ 


3.3 Equilibrium points and stability 


The next result describes the situation of the equilibrium points and their stability. Recall that the cone 
tC used to order the state space has been defined in (121 (in Theorem [ 5 ]). 


Theorem 7. System © possesses four equilibrium points, denoted £ 0 , 0 ; %u,Ch £ 0 ,w and xu,w and 
corresponding respectively to zero population, disease-free state, complete infestation, and coexistence. 
Moreover, the latter fulfill the following inequalities: 


xu ,0 xu,w <C/c £ 0 ,w and xu ,0 -C/c #o,o x o,w ■ (15) 

Last, the equilibrium points Xu ,0 and xq, w are locally asymptotically stable (LAS), while the two other 
ones are unstable. 

The proof of Theorem [T] is decomposed in the following sections. 
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3.3.1 Proof of Theorem [7]— Computation and ordering of the equilibrium points 

One here computes the equilibrium points. The latter verify 

KuTZ-o — ~v~~a — Ajj - (1 + L w + Lu)Lu = 0 (16a) 

Ay + Aw 

Iw'Tt'o' Aw — (1 + Lw + Ljj)Lw = 0 (16b) 

Ljj = 7 uAjj, L w = jwAw (16c) 


The point Xg,o := (0,0,0,0) is clearly an equilibrium. Let us look for an equilibrium xjj,o := 
(Ly, Ay, 0, 0). The quantities Ly, Ay then have to verify 


'Yu'R-o Au — (1 + Ly)Ly — 0, Ly — 7 uAy . 

Dividing by Ly ^ 0 yields 1 + Ly = TZq . One thus gets the unique solution of this form verifying 


(17) 


L'b = K -1, 


V u — 1 

a* - 172_ t 

Au ~ 

7 u 


which is positive due to hypothesis ©• 

Similarly, one now looks for an equilibrium defined as Xq : w := (0, 0, Ufa, A^)- The values of L\ v , A* w 
must verify 

Iw'R-'a A\y — (1 + Lw)Lw = 0 , Lw = lw Aw ■ 

This is identical to and as for the Xu ,0 case, one gets a unique, positive, solution, namely 


T* _ 7 pW 1 

^ W ~ /c 0 ~~ U 


A* — 

A W — 


TZ^ - 1 

lw 


(18) 


We show now that system ®> also admits a unique coexistence equilibrium with positive components 
xy,w = (Lf/*, Ay , L ^, A^). Coming back to (16) and expressing the value of the factor common to the 
first and second identity leads to 


A * * 

1 + L*j+L% = lwKj^=^ 

^w 


‘'0 


= luK[ 


A ** 

u A u 


A* * A** 

Au = nb Au 


0 A** /!** T ** 0 d** I A ** 

S±tt -r S±iy SiTT I Si 


■w 


One thus deduces 


4 ** 

A u 


K_ 

K ’ 


A * * 1 A ** 

and one can express all three remaining unknowns in function of Ayfr: 


L ** 4 ** 4 ** 

w — Iw-A-W, Ay — 


^ ... 
nU - kf lv ’ 


l* v * = luAb* = lu 


1Z W 

/V 0 A** 

rS± T, 


/v 0 

Using the value of Ly and LS now yields the relation 


- Uf W ' 


— 1 = Ly + Lw = lw ( 1 + 


lu 


7 ?w 


A** 

A w > 


lw IZq — 'Rff j 

which has a unique, positive, solution when (J5| holds. Setting for sake of simplicity 

lu 


5 := 


K 


lw TZb ~ K 


the fourth equilibrium is finally given by 

L u = - 1). 

L w = - 1 ). 


A * u * = ITTsyru^ ~ 1} 
A% = „ , * « - 1) 


(1 + 5) j w 
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We have so far exhibited all the equilibrium points. 

Notice that the last equilibrium can be expressed alternatively by use of the values of the equilibrium 
Xq , W '- 


j ** f * 

Lu - TVs Wl 

j ** _ ^ j * 


4 * * _ 


S 7 w ,* 

1 + ^ lu W 

1 


\ ** _ "*■ \ * 

Aw - rTs w 


J 5)t r * 

L/JJ V l^w/ V Ljj 


and this provides straightforward comparison result: 

and L^ < L^y < Lj 

and thus A** = 7 V L** < 7 V L* = A*, for i] € {U,W} and, therefore, 

and 


AT < a* l 






( 21 a) 

( 21 b) 

( 22 a) 

( 22 b) 


the second inequality being directly deduced from (21b). The relations (22) allow to establish the 
inequalities (15). 


3 . 3.2 Proof of Theorem [ 7 ]— Local stability analysis 

The local stability analysis is conducted through the Jacobian matrices. Recall that the gradient has 
been computed in (131. 


Stability of * 0 , 0 • The value of V/ at cco,o is not defined. However, the trajectories issued from points 
in x { 0 } 2 , resp. { 0 }~ x R+ clearly remain in the respective subspace, and the stability in these 
directions is controlled by the spectrum of the matrices 


-1 iuK\ 
1 - 7 u J 


and 


-1 7 wVT 

1 -7 w 


Their eigenvalues are 


\ (—(1 + 7r?) ± ((1 + 7r,) 2 + 47(^0 - T ) 1/2 ) » V = U , 


= u. w. 


One of each pair is positive, due to condition 0. and xgy is thus unstable. 


Stability of xjj, o- 


Using (13), the gradient Vf(xu,o) of / at xu ,0 is the 


upper block-triangular matrix 


(-1 - 2 L* v 

luV 

T* 

~ L U 

-7 uK) 



-V; 

7 uK 

l-^o 

-7 uK\ 

1 

-7 u 

0 

0 


1 

-7 u 

0 

0 

0 

0 

-1 - L* v 

7 wK 



0 

0 

—1Z U 

/c o 

iwK 

V 0 

0 

1 

-iw j 



0 

0 

1 

- 7 w J 


Using the same arguments than the ones used to study a; 0j o to assess the stability of 2 x 2 matrices, the 
two diagonal blocks are asymptotically stable if, respectively, 


- yT - < 1 and VT < V ■ 

2V -1 00 

These conditions are realized when © holds. In conclusion, the equilibrium Xjj.o is locally asymptotically 
stable. 
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Stability of xq,w The gradient Vf(xo,w) of / at Xq,w is the lower block-triangular matrix 


/-I ~ L W 

0 

0 

° ^ 


/ 

0 

0 

° \ 

1 

-7 u 

0 

0 


1 

-7 u 

0 

0 

~ L W 

0 

-1 - 2 L* w 

TwK™ 


1 - nf 

0 

1 - 2 TVS' 

7 wK™ 

V 0 

0 

1 

-7 w / 


\ 0 

0 

1 

-7 w J 


The left-upper block is a Hurwitz matrix, while the asymptotic stability of the second one is equivalent 
to 

2 - 1 > , 

that is TZq' > 1, which is true, due to hypothesis ©• The equilibrium x$ t w is thus locally asymptotically 
stable. 


Stability of xjj,w- The instability of xjj,w can be proved by showing that the determinant of the 
Jacobian matrix V f{xu,w) is negative, which, together with the fact that the state space has even 
dimension 4, establishes the existence of a positive real root to the characteristic polynomial; and thus 
that the Jacobian is not a Hurwitz matrix. This argument yields lengthy computations. 

It is more appropriate to use here the monotonicity properties of system ([6]), established in Theorem 
[5j As a matter of fact, bringing together the inequalities (15) (already proved in the end of the previous 
section, see ([22|), the asymptotical stability of Xu,o and £o,w and the strongly order-preserving property 


of the system, (Smith, 1995 Theorem 2.2) shows that the intermediary point xu,w cannot be stable. 


This finally achieves the stability analysis, as well as the proof of Theorem [7] 


3.4 Positively invariant sets and basins of attraction 

We exploit further in the following result the inequalities ordering the equilibrium points to have sup¬ 
plementary informations on some invariant sets. 

Theorem 8. The order interval 


\xu,o\ £o,w]/c '■= {a-’ € : xu ,0 </c x <k. C R 4 } 

is positively invariant for system ([6]). Moreover, the order interval \xu,w',x o,vc]ac (resp. lxu,o]Xu,w}lc) 
is contained in the basin of attraction of Xq,w (resp. Xq,u)- 


Proof. The positive invariance properties are direct consequences of the monotonicity properties exhibited 
in Theorem [5j More precisely, endowing the state space with the ordering induced by the cone K 
(see (JT 2 J) ), the autonomous system ([ 3 ]) induces a monotone flow in R 4 , strongly monotone in R+ \ 
(R^_ x {0} 2 U {0} 2 x R 2 ). As the trajectories are bounded, the set of initial points whose corresponding 
trajectories do not converge to one of the equilibria is of zero measure (Hirsch 1988 Theorem 7.8). 
Among the equilibria, only Xu ,0 and Xq,w are locally stable. 

The same rationale applies for any trajectory with initial condition in the order interval \xu,w\ ^o,vc]k:\ 
(R^_ x {0} 2 U {0} 2 x R^_), and the convergence (for almost every initial condition) in this interval can 
only occur towards x 0 ,w : the latter is therefore included in the basin of attraction. The same argument 
applies for the other equilibrium xu, o- □ 


4 Analysis of the controlled system 


4.1 A class of static output-feedback control laws 

The following feedback law will be considered in the sequel: 

u = KLjj 

for adequate (positive) values of the scalar gain K. Writing 

0 

e;= 0 

w 


(24) 
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one obtains the closed-loop system: 
or in developed form: 


x = f(x) + KBe r x , 


Lu = —^- 7 — A v - (1 + L w + Lu)Lu 

Ajj = L v - iuA v 

Lw = IwT^q Aw — (1 + Lw + Lu)Lw + A Lit 
Aw = Lw — "iwAw 


(25) 


(26a) 

(26b) 

(26c) 

(26d) 


The basic results gathered in the following theorem can be demonstrated by use of the same arguments 
than for Theorem [l] The proof presents no difficulty and is left to the reader. 


Theorem 9. For any initial value in Ri, there exists a unique solution to the initial value problem 
associated to system (25). The latter is defined on [0, +oo), depends continuously on the initial conditions 


and takes on values in Ri.. Moreover, it is uniformly ultimately bounded. 


4.2 Equilibrium points and critical gain 

We now turn to the study of the equilibrium points. The following result shows that, for gains larger than 
certain critical value, the only locally asymptotically stable equilibrium is xo,w■ Moreover, the explicit 
value of this critical number depends only upon the basic offspring numbers of the two populations and 
the ratio between their mortality rates, which are all scale-free information. 

Theorem 10. If the feedback gain K is such that 


K > IC 



(27) 


then the closed-loop system (25) possesses two equilibrium points, namely 2o,o an d x o,w> an d their local 
stability properties are not modified (i.e. Xq.q is unstable and Xo,w is locally stable). 

Proof of Theorem\10\ 

• The equilibrium points of system (25) are the points that verify 


f{x) + KBe T x = 0 . 


(28) 


Clearly, the points xg,o and Xq,w are still equilibria of system (25), as in these points e T x = Lu = 0; 
and there are no other equilibria with Lu = 0. In fact, from the third and fourth equations of (28), one 
should obtain 

0 = ^wP-^ Aw - (1 + L W )L W = ( L* w — L W )L W , 


where the second identity follows from the definition of L ^ 
equilibria than Xo : o and Xo,w- 

At any equilibrium point such that Lu ^ 0, (p8|) yields 


in (18). Let us show that there are no other 


0 = Jw'R-o Aw ~ (1 + Lw + Lu)L w + L-Lu = — 1 — Lw — Lu)Lw + PLu (29) 


and thus Lw ^ 0. At such equilibrium point, one should have 

Au 




Au + A 


= n. 


w 


w 


ts U tsW i 

A -— = A 0 + 
A w 


JfL K A u 


7 w 


A 


w 


Defining the unknown quantity 


the latter should fulfill 



n L n 


1 + 1 


= n, 


w 


K9 

lw 


(30) 
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that is 


7 u 
7 w 


K9 2 


-k + tz^-tz^ 

7 w 


e + -R™ = o 


The roots of this equation are given by 


7 w 
27 uK 


7 u 


- ( —K + 1Z% - 1Z\ 
,7 w 


±\/(^-K + n™-n\>) - 4 tv^—k 


7 w 


7 u 


7 w 


For positive values of K , there exist real nonnegative solutions to this equation if, and only if, 

2 

/ 0/7-r TTT T-r N 

and 


—K + TZ^-1Z% < 0 
7 w 


-^-K + 1Z^ - 1Z% 
7 w 


> 4B^ — K , 
Iw 


that is if and only if 

This is equivalent to 

or again 


m~K + TZ^ - nV < —2. llZ™^-K . 


7 w 


^o — 


■ 7 u 


7 w 


K < 


7 w 

iw 
7 u 


—K + J-R™ ) -vM < 0 


hz^-Jnf 1 =k* 


As hypothesis (271 is incompatible with the previous inequality, we deduce that system (25) possess only 
two equilibrium points. 

• We now study the local stability properties of the latter, by applying adequate modifications to the 
gradient exhibited in (131 and used in Section 3.3.2 to study the stability of the uncontrolled model 
equilibria. In fact, one just has to add to V/(xo,w-) the term 

/ 0 0 0 0\ 

0 0 0 0 
A' 0 0 0 
\0 0 0 0 / 

It is clear, due to the form of this additional term, that the characteristic polynomial of the system 
obtained from linearizing ( |25| ) at Xo,w is affine with respect to K, and that, for K = 0, it coincides with 
the characteristic polynomial of the linearization of ^ . 

Developing the determinant det(A/ — V/(x 0 i w) ~ KBe T ) (see (23)), the additional term is equal, for 
the feedback control law defined in (24), to 


- K 


0 0 0 

A + 7 u 0 0 

0 —1 A+ 7 w 


= 0 . 


Therefore the local behavior is not perturbed, and the asymptotic stability of the equilibrium Xo,w is 
conserved when the control term KBe T x is added. This achieves the proof of Theorem 10 □ 


Remark 11 . It can be checked from the latter proof that the previous result is not true when the effect 
of cytoplasmic incompatibility is absent. The latter is materialized by the term A — present in the 
first line in equation (|25|). When replacing this term by 1, (30) is replaced by 


V u — V 

'■Co — /■c 0 


w 


lu 
7 w 


K9 , 


which possesses the positive solution 


6 = 


7 w 
luK 


«~TZ, 


o ) 


leading to a coexistence equilibrium solution, in addition to xo,o an d %o,w- 
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4.3 Global stability issues 

We now turn to the most innovative part of this paper, namely the global behavior of the closed-loop 


system (25). The result we establish here shows that the introduction of infected larvae according to 


the proportional feedback law (241 yields conclusive infestation when the gain is larger than the critical 


value. More precisely, we have the following convergence result. 


Theorem 12. If K > K*, all trajectories of system (25) issuing from a point in distinct from Xq,o 
converge towards the complete infestation equilibrium Xq, w ■ 

Notice that strictly speaking, Theorem |12| is an almost global convergence result: convergence towards 
the complete infestation equilibrium is ensured, except for a zero measure set of initial conditions. 
However, in the present case, this set is reduced to the unstable equilibrium. 

Two attempts to prove Theorem [12] are rapidly presented in Sections |4.3.1| and |4.3.2| The main 
interest is to show how these quite natural approaches fail to provide information on the asymptotic 
behavior and, therefore, that a new approach is needed. Next a third conclusive method is exposed in 
Section [4. 3. 3[ where Theorem [12] is finally proved. 


4.3.1 Global stability of a singularly perturbed system, by LaSalle’s invariance principle 

We present now a first attempt, based on Lyapunov techniques. Consider the simpler system 

Lu)Lu (31a) 

<u -r iu^W 

Lw = 7 Lw — (1 + Lw + Ljj)Lw + KLu (31b) 


L V = K - r 1W ^ U r L U - (1 + L W 


As can be easily verified, system (31) is deduced from (26) by applying singular perturbation, formally 


putting 0 = Ljj — 'yuAjj , 0 = Lw — 7 wAw- hr other words, we assume here that (26b) and (26d| are 


fast dynamics, while (26a) and (26c) are comparatively much slower. 


The proof of well-posedness and positiveness of system (31) presents no difficulties, one states directly 
the asymptotic properties of this system. 


Theorem 13. Assume I\ > K*. Then system (31) possesses two equilibria, namely (0,0) and (0, Lw). 
The former one is unstable, while the latter one is locally asymptotically stable. Last, all trajectories of 
system (31) converges towards (0, L^), except the unstable equilibrium (0,0) itself. 


Proof. The first two assertions, proved in the same way than the similar properties were established for 
system ([ 3 ]), are not detailed. The third point is proved by considering the following candidate Lyapunov 
function: 

V(Lu,L w ) := 


Lu + Lw 

,_l_ \ {(0,0)}. Letting V(Lu , Lw) denote the value of the derivative with 


defined in the invariant set 
respect to time of V(Ljj( t), Lw(t)) along the trajectories of © , we have 


V = 


LjjLw ~ LwLjj 


(L v + Lw ) 2 


(Lu 


J W) 


r,u IwLu 
0 IwLu 


LuLw — Po LuLw + KLf 


j u 


luLw 

«LuL* w 


Now, 


(Lu + L w ) 2 (^wLu + TuLw) 

+ «Lw + Kju ~ Pq Lw) L v L w + K^ w L\f) . 


Pq' luLw + (P^lw + K'yu ~ PqIw) L v L w + K^w^nj 
\jP q TuLw - \JKlwLu 

+ ^2 yjKP^ju 7 w + P™ 7 w + P-Lu - P^Lw^j L V L W ■ 
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One verifies that the following factorization holds: 


luK + ^uIw'Sk + 7 w - TZ^lw 

- (v^V^))(^ + )f 

= lu (yK-Vi^)(VK + ^(yi^+y/K^ 



(32) 


where K* is the critical gain given in (271. Obviously the quantity in the right-hand side of (321 is 
positive whenever I\ > K*. One thus has 


V{Lu, Lw) < 0 


for any ( Ljr,Lw ) G R+ \ {(0,0)}, with equality if and only if Ljj = 0. LaSalle’s Invariance Principle 
(La Salle| 19T6j) is then used to conclude. □ 


While the Lyapunov function L J^ r L used in the previous proof is quite appealing, it has not proved 
possible to extend this idea to the complete controlled system (25). 


4.3.2 A monotone control system perspective 

The second method now explored is an attempt to apply the results on monotone control systems , as 
worked out in particular by |Angeli and Sontag (2003), see also the works by Gouze (1988); Cosner ( 1997| ) ; 
Enciso and Sontag (2006); {Enciso (2014). The principle of this approach consists in decomposing the 
system under study as a monotone input-output system with feedback. Given the fact (see Theorem [5]) 
that the uncontrolled system is monotone, a most natural way to do this is to write system (25) as 


x = f(x ) + KBu, y = Ljj = e T x 
u = y 


(33a) 

(33b) 


Arguing as in Theorem[5j one can establish that the input-to-state map u4i given by (33a) is monotone 


when the state space is endowed with the ordering <x\ while the state-to-output map x t->- y = e T x = Ljj 
is anti-monotone. We are thus in the configuration of a so-called monotone system with negative feedback. 
In such a case, the study of asymptotics of the system obtained when closing the loop by the unitary 


feedback (33b) can be done by introducing static characteristics (Angeli and Sontag 2003 Enciso 2014). 


By definition, when it exists, the input-state characteristic kx associates to any constant input u the 
corresponding value kx{u) of the unique globally asymptotically stable equilibrium; and the input- 
output characteristic k, obtained by composing kx with the state-to-output map, associates to u the 
corresponding output value k(u). Notice that for a monotone system with negative feedback (as it 
is the case here), the map k is non-increasing. The inspiring results demonstrated in the references 
above establish general conditions under which existence and stability of fixed points of the input-output 
characteristic k permit to deduce the convergence of every (or almost every) trajectory of the closed-loop 
system towards a locally asymptotically stable equilibrium; and that these equilibria are in one-to-one 
correspondence with the stable fixed points (through the state-to-output map). 

However, things become immediately complicated in the case of system ( |33[ ): constant input u = 0 
leads to the uncontrolled system (§, which has been proved to possess two locally asymptotically stable 
equilibria (and two unstable ones). In such a case, one can consider multivalued input-state and input- 
output characteristics, as made by Malisoff and Leenheer| ( |2006[ ), or |Gedeon a nd Hines| (J2009). The 
corresponding input-output characteristic is drawn in red in Figure [l] (for the parameter values given in 
Section [5] below) and for the corresponding critical value K* of K. 

As seen in the figure, the input-output characteristic has basically two branches (both drawn in red). 
The first one merges with the horizontal axis: it corresponds to a branch of equilibria with null value of 
the output y = Ljj that departs from xo : w for u = 0. The second one is a decreasing curve, defined for 
values of u ranging from zero to a value close to 2.32: it corresponds to the output value of a branch of 
equilibria departing from xjj.o- The blue curve, which does not pertain to the input-output characteristic, 
indicates the output values of a branch of unstable equilibria originating from xjj : w and that vanishes 
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Figure 1: 
(33a I (in 


The multivalued input-output characteristic u >->• k(u) corresponding to input-output system 
red) and the diagonal u K > u (in black), in the critical case K* = K. See text for explanations 


together with the upper curve. The diagonal line (that determines the fixed points of k) is also shown. 
It is tangent to the blue curve, due to the fact that K = K* here. For smaller values of K , the diagonal 
intersects twice the two upper branches; while for larger values of K the only intersection between the 
input-output characteristic and the diagonal is the origin. 

For K > K *, the complete infestation equilibrium is therefore the only fixed point of the multivalued 
map k. But the iterative sequences Uk+i = fc(ufe) do not converge systematically towards this point. In 
fact, the only information that can be deduced from the results applicable to cases of multivalued input- 
output characteristic (Malisoff and Leenheer 2006| Gedeon and Hines 20091, is that all trajectories are 
bounded, and that the output y = Ljj fulfills the following inequalities: 


0 < limsupy(t) < A; sup (liminfy(t) 

V t— H-oo 


t —^ -|-oo 


Here fc sup denotes the discontinuous function whose curve is equal to the upper branch until it vanishes, 
and merges afterwards with the horizontal axis. 


As a conclusion, the decomposition (33), that seemed a natural framework to analyze the behavior 


of the controlled system (25) immediately fails to produce a global vision of the asymptotic behavior. 


4.3.3 Monotonicity revisited and proof of the global stability 

We now concentrate without supplementary detour on the proof of Theorem 


in working on an alternative decomposition of system (251, different from (33). Define first 


12 The principle consists 



f z if z < 0 

f z if z > 0 


\z _ := < and 

|*|+ := 


1 0 otherwise 

1 0 otherwise 

Clearly, one has 

0 = \z\- + M + , 

zeK. 


(34) 
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We introduce the following decomposition that will show convenient in the proof of Theorem 12 


Lu = 


Au 


0 A u + A 


-Au — (1 + Lw + Lu)Lu 


U T -fi w 
Au = L v - 'YuAu 

Lw — IwLto l Aw — (1 + Lw)Lw + |K — Lw\-Lu + Ku 
Aw = Lw — IwAw 
Lw 


V = 


1 - 


K 


Lu 


(35a) 

(35b) 

(35c) 

(35d) 

(35e) 


u = y indeed yields system (251. 


As a matter of fact, using property (34), one sees easily that the closing of the input-output link (35) by 


Next, we state and prove Lemmas 14 15 and 16 that will be used in the proof of the main result, 
Theorem fl2l 


Lemma 14. 


(35). 


For any integrable u taking on nonnegative values, the set is positively invariant by 


Proof of Lemma [77 
border of 


The key point is that \K — Lw\-Lu = 0 when 0 < Lw < K. Therefore, near the 
\ where Lw = 0, the system (35) behaves locally as Lw = ^wKq' Aw — (1 + Lw)Lw + Ku. 
The fact that Lw > 0 whenever Lw = 0 then forbids escape from the set R% by this side. The same 
happens for the other three variables: their derivatives are nonnegative at the points where they vanish. 
Hence, the trajectories can neither escape by the other sides. This establishes the positive invariance of 
R\ and achieves the proof of Lemma [l4[ □ 


Lemma 15. The input-output system (35) is monotone with negative feedback, when the state space is 


endowed with the order defined in Theorem [5[ 


Proof of Lemma 15. The right-hand sides of (j35aj , (35b I and (35d) have been studied in Theorem [5j 
The right-hand side of (35c) is clearly increasing with respect to Aw and u, and non-increasing with 


respect to Lu- So the input-to-state map is monotone when the state space is endowed with the order 
defined in Theorem [5] 


On the other hand, the state-to-output map defined by (35e) is non-increasing with respect to Lw, 


and non-decreasing with respect to Lu- Therefore, it is anti-monotone with respect to the ordering used 
in the state space. This achieves the proof of Lemma [15} □ 


Lemma 16. For any constant nonnegative input, system (35a) —(35d I possesses a unique LAS equilib¬ 
rium. The latter yields null value of Lu- 


Moreover, the solution of the input-output system (35a) —(35d) converges towards the corresponding 


equilibrium when time qoes to infinity (and in particular the output Ltt converqes to zero), except possibly 
ifu = 0 and L w (0) = 0, A w (0) = 0. 


Lemma 16 does not allow to define in the usual way an (identically null) input-output characteristic: 
when u > 0, the solution of the input-output system (35) converges to the LAS equilibrium for any initial 
condition, and is such that 

lim Lu{t ) = 0 ; 

t —>-+oo 

but when u = 0, this property is only guaranteed if (Lw(0 ), 4^(0)) ^ (0, 0). 

However, Lemma [16] allows to define a weaker notion of characteristic called input-output quasi-cha¬ 


racteristic (Angeli et al 


2004). Contrary to the stronger notion (presented in Section 4.3.2), for any 


constant input value, the convergence to the asymptotic value does not have to be global: convergence 
is only required to occur in a full measure set (i.e. a set whose complementary is of zero-measure). This 
is exactly what happens in the present situation, as the set of those x G R% for which Lw{0 ) = 0 and 
Aw(0 ) = 0 is negligible. 


As a consequence, one can in fact show, with the tools developed by Angeli et al (2004), that Xq,w 


is a/raost-globally attractive for system (25), i.e. that it attracts all trajectories, except possibly those 


departing from certain zero-measure set. The difficulty we now face is that this negligible set contains 
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a priori all the points ofR% such that (Lw{0), Aw (0)) = (0,0), and in particular the equilibrium xjj, o- 
However, it is quite natural to assume that, at the moment where the release begins, the system departs 
from a Wolbachia- free situation. In order to ensure that these situations too are concerned by the 
convergence to Xo t w> more precise arguments are therefore needed. 

We first prove Lemma |16[ and demonstrate in the sequel that, as announced in Theorem |12[ the 
convergence occurs for all trajectories, except the unique trajectory immobile at Xg,o- 


Proof of Lemma \Tf\ One first studies the equilibria of system (35a) (35d), for constant inputs u : 
u(t) = u , for u £ R + . Clearly, the set of these equilibria is the union of two sets: the set of equilibria of 
(35aJ + (35b)+(35d) and 

Lw = "IwLZ^ Aw ~ (1 + Lw)Lw + Ku (36a) 


such that I\ — Lw > 0; and the set of equilibria of (35a)+(35b)+(35d) and 

Lw = Aw — (1 + Lw + Ljj)Lw + KLu Ku 


(36b) 


such that K — Lw < 0. 


Consider first the system (35al + (35bl+(36al+(35d|. As can be seen, the latter is decoupled, since 


Ljj is not anymore present in the right-hand side of (36a). One shows without difficulty that there exists 
a unique equilibrium in Kl, which is LAS and characterized by 


L w — 


Tl™ - 1 + - l) 2 + 41<uj 


(37) 


provided that this expression verifies Lw < K. Another equilibrium exists, which is Xq.o if u = 0, but 
which has negative value of Lw if u > 0, and is therefore discarded, due to Lemma [I4| 


Consider now the second case, of system (35a)+(35b)+(36b)-|-(35d). Arguing as in the proof of 


Theorem 10 the only equilibria that may exist are such that Ljj = 0. As a matter of fact, for a solution 


with nonzero Lu , a term Ku in the right-hand side of (29) could be written, jointly with KLjj, as K'Lu 
for some K' > K > K*, leading therefore to Ljj = 0 and a contradiction. Therefore, any potential 
equilibrium has to fulfill Ljj = 0, and the only possibility is given by (37) if this expression verifies 
Lw > K. 

Putting together the two cases, one sees that: 

* there exist two equilibria, Xq,o and xg, w, if u = 0; 

* there exists a unique equilibrium if u > 0. 

* In any case, the corresponding output value is 0. 

Now, for any constant input u{t) = u, system (35) is strongly order-preserving, just as system © 
was shown to be (Theorem [5]). Then, the uniqueness of equilibrium in the case where u > 0 allows to 
Theorem 2.3.1, p. 18) and to deduce that all trajectories in converge to this unique 


use (Smith 


equr 


1995 


ibrium. When u = 0, applying (Smith 1995 Theorem 2.2.1, p. 17) shows that every trajectory 


converges towards one of the two equilibria £ 0,0 and Xq,w- The behavior of system (351 in the vicinity 
of xo t o obeys the equations 


Lw = IwKq Aw — (1 + Lw)Lw, Aw = Lw — "twAw ■ 

This system is monotone and the projection of Xq,w attracts all trajectories, except if Lw( 0) = 0 and 
Aw{ 0) = 0. This achieves the proof of Lemma 16 □ 


One is now ready to achieve the proof of Theorem [l2| We define y(t-,x 0 ,u) the output of system 


(35a) (35d) corresponding to the input signal u and the initial state value x$. For any trajectory of the 


closed-loop system (25), we will denote indifferently u and y, in order to exploit the formalism of the 


input-output decomposition given in (35). 

First of all, recall that, due to Theorem |9j all trajectories of (25) are bounded. Therefore, for any 
nonnegative initial condition Xg, 


0 < lim inf y(t; xg,u) < lim sup y(t] xg, u ) < +00 Wxg £ 

t-H-oo t-H-00 


| 4 + 


(38) 
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Assume first 


liminf y(t; Xq, u) > 0 . 
£—>■+00 


(39) 


Using monotonicity of the input-output system to compare trajectories with different inputs, one obtains 
from the fact that 


that 


Using Lemma 


Ve > 0,3 T >0, t > T => u(t) > lim inf u(t') — e 

£'—>■ +00 


Ve > 0, lim sup y(t; Xq , it) < lim sup y (t; xq, lim inf u(t') — e 


£—>■+00 


£—^+00 


>+00 
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for e > 0 smaller than liminf y(t'; Xq,u) yields 

t'—H- 00 


lim sup y ( t;x 0 , lim inf u(t') — £ ) = 0 . 
£—>■+00 \ t'^+oo J 


By putting together (38), (401 and (41), one gets: 


0 < liminf y(t; Xo, u) < limsupy(t; Xg, u) < 0 , 

t-H-00 t-y+o o 


which is absurd. This shows consequently that the premise (39) was erroneous. 
We thus have 


for all trajectories. Assume 
As above, one can deduce that 


liminf y(t\ Xq, u) =0 
£—>+00 


(L w (0),A w (0)) / (0,0) . 


(40) 

(41) 

(42) 

(43) 

(44) 


£—>•+00 

and therefore that 


0 < liminf y(t\ xq, u) < \inisupy{t',Xo,u) < lim sup y t; Xq, liminf u(t') 1=0, 

x . , x , . \ £'_>■ + OO ' 


£—>•+00 


£—>+00 


lim y[t ; xq, u) = 0 . 
£—>•+00 


Now Lemma 16 permits to deduce from (45) and (44) that 

lim x(t; Xq,u) = Xo w ■ 
£—>+00 


(45) 


(46) 


On the other hand, if (44) is false but (Lj/(0),A+O)) 7 ^ (0,0), then it is easy to show that 
(Lw(t), Aw(t)) 7 ^ (0,0) for some t > 0 (and indeed for any t > 0). As a matter of fact, due to the 
presence of the control term (which is continuous and initially positive), L\y is certainly positive on a 
sufficient small punctured open neighborhood of t = 0. This in turn yields the same property for Aw, 
due to the linearity of its evolution. The analysis previously conducted in the case where (44) is true, 
can therefore be applied in the present case (where (44) is false but (Ljj( 0), A+O)) 7 ^ (0, 0)) from a new, 
positive, initial time instant. It allows to conclude similarly that (45) and (46) hold. 

As a conclusion, the convergence to Xq,w occurs in any case, except if (Lvv(0), Aw (0)) = (Ljj(0), Ay ( 0 )) 
= (0,0), that is except if a;(0) = £ 0 , 0 • This achieves the proof of Theorem 
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5 Numerical simulations 

We present some illustrative simulations, with the following realistic values: 

lu = 0-8, 7 w = 1, K = 5, 7C = 3 

Notice that the mortality is higher for the Wolbachia infected population (yj/ < ^w), and its sustain¬ 
ability is inferior (7+ > 7?+). The critical gain value can be computed and is equal to 

I\* ~ 0.318 

Figures [ 2 ] to [5] show the evolution of the state variables L v (in blue) and L w (in en) as functions 
of time. The initial state is the Wolbachia-iree equilibrium xu,o, and the gain values are respectively 
chosen to be 1, 0.5, 0.35 and 0.3. The last value, smaller than the critical value K *, yields convergence 
to a coexistence equilibrium. 
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Figure 2: Evolution of Ljj(t) and L\y(t) as functions of time, for K = 1 



Figure 3: Evolution of Ljj(t) and as functions of time, for K = 1 and K = 0.5 


6 Conclusions and further studies 

We presented and analyzed a model for the infestation by bacterium Wolbachia of a population of 
mosquitoes — typically one of the genera Aedes involved in the transmission of arboviroses such as 
yellow fever, dengue fever or chikungunya. A method of implementation based on the introduction of 
a quantity of contaminated insects proportional to the size of the healthy population was proposed and 
shown, analytically and by simulation, to be capable to spread successfully the bacteria provided the 
gain is sufficiently large. This feedback method requires continuous measurement of the population. Its 
main interest with respect to the release(s) of a predefined quantity, is the reduction of the number of 
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Figure 4: Evolution of Ljjit) and Lw(t) as functions of time, for K = 1, K = 0.5 and K = 0.35 



Figure 5: Evolution of Ljj(t) and Ly\ (t) as functions of time, for K = 1, K = 0.5, K = 0.35 and 
K = 0.3. The last value is smaller than the critical value K* ~ 0.318, and a coexistence equilibrium 
appears asymptotically 


released mosquitoes, and thus of the treatment cost, without jeopardizing the success of the introduction 
of the bacteria — something which can happen e.g. in case of underestimation of the initial population 
size. To our knowledge, this is the first use of the control theory notion of feedback in such a context. 

Among other steps leading to application, the adaptation to effective conditions has to be done. First, 
the model presented here has been chosen continuous in time for simplicity, but passing to discrete-time 
system seems to present a priori no difficulties. Also, the present framework assumes measurement of a 
larva stage of the healthy portion of the population, and as well release of larva stage of the contaminated 
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one. The practical conditions can be different, and the method can be adapted in consequence (leading 
though to similar, but different, convergence questions). Last, issues of robustness with respect to the 
uncertainties of the parameters that describes the dynamics have not been tackled here. 

An advantage of the present modeling framework is to open the way to comparisons with optimal 
policies — for example the one that minimizes the total number of released mosquitoes, while succeeding 
in spreading Wolbachia. This point will be studied in a next future. Also, this framework provides a 
first basis to consider questions related to strategy improvement by mitigating several control principles, 
or to the complex phenomena of interaction between different vector species and different arboviruses, 
that may occur in the context of control of different diseases. 

From a mathematical point of view, one of the difficulties of the study is that the system presents 
two stable equilibria, corresponding to Wolbachia-free situation and complete infestation. While the 
key arguments are based on the theory of input-output monotone systems developed after Angeli and 


Sontag (2003), none of the posterior refinements to multivalued characteristics or quasi-characteristics 


allowed to establish formally the main convergence result, and adequate adaptation had to be achieved. 
Extensions in this direction are presently studied. 
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